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Final  Report 
Tai  C.  Chiang 
Abstract 

This  proposal  is  a  study  of  the  electronic  structure  and  quantum  effects  of  layer  systems  based  on 
hexagonal  silicon  carbide,  a  large-gap  electronic  substrate  material  suitable  for  high-power  and 
radiation-resistant  applications.  High  temperature  processing  of  SiC  surfaces  in  vacuum  leads  to 
reconstructed  surfaces  with  varying  chemical  compositions  and  atomic  organizations.  Prolonged 
heating  depletes  Si  from  the  surface,  resulting  in  the  formation  of  a  graphene  sheet  or  a  graphite¬ 
like  structure  made  of  multiple  graphene  sheets.  These  systems  possess  useful  and/or  unusual 
electronic  structures  including  Dirac  fermion  dispersions  near  the  Fermi  level  for  a  graphene 
sheet,  band  splittings  in  multiple  sheets,  and  band  bending  within  the  SiC.  The  potential  for 
novel  electronic  applications  has  driven  large  prevailing  research  efforts  within  the  materials 
science  community.  Our  research  has  focused  on  these  low-dimensional  systems,  the  formation 
of  ultrathin  metallic  overlayers,  and  the  resulting  quantum  physical  properties  of  the  composite 
systems.  Key  discoveries  and  advances  include: 

Mapping  by  angle-resolved  photoemission  spectroscopy  of  the  spectral  functions  of  graphite  and 
graphene  layers  at  low  temperatures  reveals  a  heretofore  unreported  gap  of  ~67  meV  at  normal 
emission.  This  gap  persists  to  room  temperature  and  beyond,  and  diminishes  for  increasing 
emission  angles.  We  show  that  this  gap  arises  from  electronic  coupling  to  out-of-plane 
vibrational  modes  at  the  K  point  in  the  surface  Brillouin  zone  in  accordance  with  conservation 
laws  and  selection  rules  governed  by  quantum  mechanics.  Our  study  suggests  a  new  approach 
for  characterizing  phonons  and  electron-phonon  coupling  in  solids. 

Electronic  chirality  near  the  Dirac  point  is  a  key  property  of  graphene  systems.  The  direction  of 
the  isospin  vector  in  monolayer  graphene  tracks  the  rotation  of  the  electronic  crystal  momentum, 
resulting  in  a  Berry's  phase  of  180  degrees.  By  contrast,  the  isospin  in  bilayer  graphene  rotates 
twice  as  fast,  resulting  in  a  Berry's  phase  of  360  degrees.  We  show  that  these  chiral  relationships 
are  connected  to  the  spectral  intensity  patterns  as  measured  by  angle-resolved  photoemission 
spectroscopy  under  different  polarization  conditions.  Specifically,  the  semicircular  patterns  at  30 
eV  for  monolayer  (bilayer)  graphene  rotate  by  +/- 90°  (+/-45°)  in  changing  from  linearly  to 
circularly  polarized  light.  Calculations  confirm  that  these  effects  arise  from  the  interplay  between 
electronic  chirality  and  coherent  interference  of  dipole  transitions  from  the  s  and  p  polarization 
components.  The  measured  circular  dichroism  thus  yields  definitive  information  about  the  phases 
of  the  electronic  states. 

Work  has  also  been  carried  out  to  map  the  phonon  dispersion  relations  of  graphene  systems  over 
wide  regions  in  the  surface  Brillouin  zone.  The  quantum  electronic  structure  of  metal  overlayers 
is  currently  under  investigation. 
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Phonon-Induced  Gaps  in  Graphene  and  Graphite  Observed  by  Angle-Resolved  Photoemission 
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Mapping  by  angle-resolved  photoemission  spectroscopy  of  the  spectral  functions  of  graphite  and 
graphene  layers  at  low  temperatures  reveals  a  heretofore  unreported  gap  of  —67  meV  at  normal  emission. 

This  gap  persists  to  room  temperature  and  beyond,  and  diminishes  for  increasing  emission  angles.  We  show 
that  this  gap  arises  from  electronic  coupling  to  out-of-plane  vibrational  modes  at  the  K  point  in  the  surface 
Brillouin  zone  in  accordance  with  conservation  laws  and  selection  rules  governed  by  quantum  mechanics. 

Our  study  suggests  a  new  approach  for  characterizing  phonons  and  electron-phonon  coupling  in  solids. 


DOI:  1 0. 1 1 03/Phy  sRevLett.  1 05. 1 36804 

Angle-resolved  photoemission  spectroscopy  (ARPES) 
is  a  powerful  technique  for  determining  the  electronic 
structure  of  materials,  and  it  has  been  applied  to  numerous 
systems  of  scientific  and  technological  interest  [1].  The 
measured  spectral  functions  are  related  to  the  quasiparticle 
self-energies,  which  are  influenced  by  interactions  with 
elementary  excitations  including  phonons.  Electron- 
phonon  coupling  is  of  special  interest  as  it  plays  a  key 
role  in  the  formation  of  charge  density  waves  and  super¬ 
conducting  states.  Its  main  signature  is  the  “kinks”  in 
dispersion  relations  near  the  Fermi  level,  which  have 
been  well  documented  for  a  variety  of  systems  including 
graphene  layers  of  various  thicknesses  and  configurations 
[2-6].  While  strong  interest  in  the  electronic  structure  of 
these  graphitic  materials  has  driven  extensive  ARPES 
studies  [2,3,7-10],  prior  work  has  mostly  focused  on  the 
quasiparticle  band  dispersion  relations  associated  with  the 
Dirac  cones.  Largely  unexplored  are  spectral  regions  far 
away  from  the  quasiparticle  bands,  where  direct  emission 
from  the  quasiparticles  is  forbidden,  but  indirect  emission 
through  coupling  to  phonons  is  allowed.  Our  work  shows 
that  ARPES  measurements  within  the  forbidden  region 
can  yield  valuable  information  about  high-order  processes 
such  as  electron-phonon  coupling.  Specifically,  the  experi¬ 
ment  reveals  a  gap  of  —67  meV  at  low  temperatures  for 
several  graphitic  materials  including  highly  oriented  pyro- 
litic  graphite  (HOPG),  natural  single-crystal  graphite 
(NSCG),  and  few-layer  graphene  (FLG)  systems  prepared 
on  SiC(0001).  A  detailed  study  of  the  temperature  and 
angular  dependences  of  the  gap  suggests  that  the  gap  arises 
from  out-of-plane  vibrational  modes  at  the  K  point,  which 
dominate  the  coupling  to  electrons  at  normal  emission. 

Our  ARPES  measurements  were  performed  at  the 
Synchrotron  Radiation  Center,  University  of  Wisconsin- 
Madison,  using  the  PGM- A,  U-NIM,  and  Wadsworth  beam 
lines.  The  energy  and  momentum  resolutions  were  20  meV 
and  0.01  A- 1 ,  respectively.  All  spectra  were  taken  with  the 
sample  temperature  at  60  K,  unless  otherwise  stated.  The 
results  from  these  three  samples  are  similar;  this  is  not 


PACS  numbers:  73.22.Pr,  71.38.-k,  79.60.Jv 

surprising  since  graphite  is  essentially  a  stack  of  graphene 
layers  weakly  bonded  together  by  van  der  Waals  forces. 

Figure  1  shows  the  electronic  band  structures  and 
phonon  dispersion  relations  for  graphite  and  graphene 
obtained  from  first-principles  calculations  based  on  the 
ABINIT  program  [11-14].  These  calculations  based  on  the 
LDA  approximation  yield  a  fairly  good  overall  description 
of  the  experiment;  more  accurate  results  based  on  the  GW 
approximation  are  available  in  the  literature  [15,16]. 
The  calculated  results  are  very  similar  for  graphene  and 
graphite  within  the  basal  plane.  The  dispersions  for 
graphite  are  small  along  the  z  direction  (perpendicular  to 


(a)  Graphite  (b)  Graphene 


(c)  Graphite  (d)  Graphene 


FIG.  1  (color  online).  Band  structures  and  phonon  dispersion 
relations  of  graphite  (a),  (c)  and  graphene  (b),  (d)  from  first- 
principles  calculations.  The  inset  in  (a)  shows  schematically  the 
two-dimensional  Fermi  surface  of  graphite  and  graphene,  which 
consists  of  6  Dirac  points  at  K  and  K' . 
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the  basal  plane).  Of  special  interest  are  the  phonons  at  the 
K  point  of  graphene  and  along  K-H  in  graphite.  Their 
energies  are  essentially  identical.  For  simplicity,  the  dis¬ 
cussion  below  will  mostly  refer  to  the  graphene  dispersion 
relations.  Generalization  to  graphite  is  straightforward. 

The  ARPES  spectrum  at  normal  emission  from  HOPG 
at  60  K  [Fig.  2(a)]  shows  a  peak  derived  from  the  i t  band 
at  about  —8.5  eV.  Emission  at  angles  near  the  K  point 
[Fig.  2(b)]  shows  a  quasi-Dirac  cone  from  the  tt  band.  The 
random  in-plane  orientation  of  graphite  crystallites  in 
HOPG  dilutes  the  emission  intensity  from  the  Dirac  cone, 
making  this  feature  very  weak.  Detailed  spectra  at  normal 
emission  and  near  K  [Fig.  2(c)]  reveal  a  gap  near  the  Fermi 
level  at  normal  emission,  as  evidenced  by  a  comparison  to 
the  spectrum  from  a  reference  Ag  film.  The  spectrum  at 
normal  emission  is  very  weak  as  indicated  by  the  relative 
scales.  It  actually  contains  a  tiny  Fermi  edge,  which  is 
likely  caused  by  defects  [17].  The  inset  shows  the  spectrum 
at  normal  emission  symmetrized  about  the  Fermi  level  (see 
supplementary  document  [18]),  which  clearly  shows  a  gap 
of  about  ±67  meV  in  an  approximately  linear  spectrum. 
This  gap,  independent  of  the  photon  energy,  becomes  less 
distinct  at  higher  temperatures  [Fig.  2(d)]  but  remains 
visible  at  room  temperature.  Similar  results  are  obtained 
from  a  NSCG  [Fig.  2(e)],  but  there  is  a  somewhat  higher 
background  in  the  gap  likely  caused  by  impurities  and 
defects  in  natural  crystals. 


FIG.  2  (color  online).  ARPES  results  from  HOPG  [(a),  (b),  (c), 
(d)]  and  NSCG  (e).  (a)  Normal  emission  spectrum  of  HOPG. 
(b)  Measured  dispersion  of  HOPG  near  K  along  a  radial  direc¬ 
tion  pointed  away  from  T.  (c)  Spectra  at  various  in-plane  mo¬ 
menta  k  for  HOPG  and  a  reference  spectrum  from  Ag.  The  inset 
is  a  zoom-in  view  of  the  gap  near  the  Fermi  level,  obtained  by 
symmetrizing  (labeled  S)  the  normal  emission  spectrum  with 
respect  to  the  Fermi  level,  (d),  (e)  Comparison  between  fits  (red 
[medium  gray]  curves)  and  normal  emission  spectra  for  HOPG 
(d)  and  NSCG  (e).  The  vertical  dash  lines  indicate  positions  of 
the  Fermi  level  and  the  edge  of  the  gap.  The  zero  level  for  each 
spectrum  is  indicated  by  a  tick  mark  on  the  right  vertical  axis. 
The  photon  energy  used  was  22  eV.  The  sample  was  at  60  K  in 
all  cases  except  for  those  labeled  otherwise. 


A  similar  gap  is  also  observed  in  FLG.  The  number  of 
monolayers  (MLs)  of  graphene  can  be  deduced  from  the 
Dirac-cone  feature  near  the  K  point  [Fig.  3(a)]  [8,9]. 
Emission  from  the  Si  2 p  core  level  [Fig.  3(b)]  decays 
rapidly  for  increasing  graphene  coverages  because  of  a 
short  photoelectron  escape  depth  [19].  At  normal  emission, 
the  6V3  X  6%/3  surface  exhibits  a  broad  surface  state  at 
about  —  0.6  eV  [Fig.  3(c)]  [20].  It  becomes  an  interface 
state  under  graphene  coverage,  and  its  intensity  is  attenu¬ 
ated  rapidly  with  increasing  graphene  coverages.  At  small 
graphene  coverages,  the  long  tail  of  the  interface  state 
enters  into  the  gap  region  and  interferes  with  line  shape 
analysis.  As  a  result,  the  gap  is  not  obvious  at  1  ML,  but  it  is 
quite  clear  for  the  thicker  graphene  layers  (2-4  MLs).  In 
each  case,  there  is  a  metallic  Fermi  edge,  likely  due  to 
defects.  The  gap  is  clearly  seen  in  the  enlarged  views  of 
the  3  and  4  ML  data  and  their  symmetrized  versions 
[Fig.  3(d)].  Also  included  is  a  spectrum  for  the  3  ML 
sample  near  the  K  point,  which  shows  no  gap.  The  gap 
feature  for  the  4  ML  case  is  somewhat  weaker,  possibly 


(b) 


(c) 


Energy  (eV) 


k-K(k') 


Energy  (eV) 


FIG.  3  (color  online).  ARPES  result  from  FLG  on  6H-SiC. 
(a)  Evolution  of  the  Dirac-cone  features  near  the  K  point  mea¬ 
sured  along  a  direction  perpendicular  to  T  K .  (b)  Si  2p  core  level, 
(c)  Normal  emission  spectra  as  a  function  of  graphene  layer 
thickness,  (d)  Normal  emission  spectra  and  symmetrized  versions 
(labeled  S)  near  the  Fermi  level  for  3  and  4  ML  graphene  layers. 
Also  shown  for  comparison  is  a  spectrum  near  the  K  point.  The 
zero  level  for  each  spectrum  is  indicated  by  a  tick  mark  on  the 
right  vertical  axis.  The  photon  energy  used  was  50  eV  in  all  cases, 
except  (b)  for  which  130  eV  photons  were  used. 


136804-2 


PRL  105,  136804  (2010) 


PHYSICAL  REVIEW  LETTERS 


week  ending 
24  SEPTEMBER  2010 


caused  by  increased  roughness  at  higher  graphene  cover¬ 
ages  (see  [18])  [7]. 

Because  of  the  large  band  gap  around  the  Fermi  level  at 
the  zone  center  (Fig.  1),  the  observed  gap  feature  must 
involve  higher  order  processes.  The  only  states  nearby  in 
energy  are  near  the  K  point.  Momentum  conservation  for 
normal  emission  from  these  states  requires  the  involvement 
of  either  an  elementary  excitation  near  the  K  point  or  a 
defect.  Defect  scattering,  being  strictly  energy  conserving, 
cannot  account  for  the  gap.  Phonon  scattering  is  the  only 
plausible  explanation. 

The  initial  state  at  the  K  point  for  graphene  is 

*Ak(f)  a  -  R)  exp(/K  •  R),  (1) 

R 

where  </>(r)  is  a  carbon  pz  orbital,  K  is  the  wave  vector  at 
the  K  point,  and  the  summation  is  over  all  lattice  vectors  R. 
The  orbital  </>( r)  is  of  even  parity  within  the  basal  plane. 
There  are  actually  two  degenerate  states;  the  other  one 
involves  a  pz  orbital  centered  about  the  other  atom  within 
the  unit  cell.  Photoexcitation  of  i/fK(r)  leads  to 

<Ak(r)  *  2>'(r  -  R)  exp(iK  •  R),  (2) 

R 

where  <$>'  is  an  excited  Wannier  function.  The  excitation 
involves  the  usual  dipole  transition  term  and  a  “surface 
transition”  term  arising  from  a  nonzero  V  •  A,  where  A  is 
the  vector  potential  [21,22].  Both  terms  are  important  for 
three-dimensional  systems,  but  only  the  dipole  term 
preserves  the  crystal  momentum  along  z.  Graphite  is  a 
quasi-two-dimensional  crystal  with  a  large  perpendicular 
lattice  constant  c  =  6.7  A  compared  to  the  z  extent  of  the 
pz  orbital  of  about  a  =  1.4  A;  the  dipole  term  is  scaled 
down  by  a  factor  of  — (a/c)2.  As  a  result,  surface  transi¬ 
tions  dominate.  Since  this  involves  a  z-dependent  interac¬ 
tion,  the  in-plane  parity  of  the  state  is  preserved.  Thus,  <fi' 
is  also  of  even  parity  within  the  basal  plane. 

The  state  if/fK  can  couple  via  a  phonon  of  wave  vector 
—  K  to  a  state  at  the  zone  center: 

*Ao (r)  «  5>"(r  -  R),  (3) 

R 

where  cf)"  is  another  excited  Wannier  function.  This  state 
can  couple  to  a  plane  wave  in  vacuum,  leading  to  normal 
emission,  if  <j>"  is  of  even  parity  within  the  basal  plane  [23]. 
Within  the  usual  rigid-ion  approximation,  the  electron- 
phonon  coupling  is 

A H  £>(0"(r)|e  •  W(r)|0'(r)),  (4) 

where  V  is  the  ionic  potential,  e  is  the  phonon  polarization 
vector,  and  Q  is  the  phonon  normal  mode  coordinate. 
Because  both  (£>’  and  4>"  are  of  even  parity  within  the  basal 
plane,  only  phonons  polarized  along  z  can  contribute  to  this 
process. 

Four  phonon  modes  of  graphene  at  the  K  point  with 
energies  clustered  around  150  meV  have  in-plane  polar¬ 


izations  (see  Fig.  1)  [24,25].  The  remaining  two  degenerate 
phonon  modes  at  energy  E 0  =  67  meV  are  polarized  along 
z,  and  these  are  the  relevant  ones.  At  low  temperatures, 
phonon  absorption  is  quenched.  Photoemission  assisted 
by  phonon  emission  should  exhibit  a  gap  of  67  meV  below 
the  Fermi  level  EF  by  energy  conservation.  At  finite  tem¬ 
peratures,  the  photoemission  intensity  becomes 

1(E)  a  (A  —  BE)^F(E  -  Ef  +  E0;  T) 

+  exP (j^)nE  ~ef-e0;T) ],  (5) 

where  F  is  the  Fermi-Dirac  distribution  function,  and  the 
factor  exp(— E0/kBT)  is  the  branching  ratio  between 
phonon  absorption  and  emission  processes.  The  effective 
(phonon-assisted)  density  of  states  and  cross  section 
variations  are  assumed  to  be  linear  (with  A  and  B  being 
constants).  At  high  temperatures,  phonon  absorption  pro¬ 
cesses  could  fill  in  the  gap.  However,  with  E0/kB  = 
770  K,  the  gap  remains  visible  even  at  room  temperature. 
The  red  (medium  gray)  curves  [Fig.  2(d)]  are  fits  to  the  data 
at  temperatures  T  =  60  and  300  K  using  Eq.  (5)  supple¬ 
mented  by  a  metallic  Fermi  edge  to  simulate  defect  scat¬ 
tering.  The  fits  yield  E0  =  64  meV,  in  good  agreement 
with  the  known  value  of  67  meV.  Analysis  of  the  other 
cases  [e.g.,  Fig.  2(e)]  yields  similar  values  of  E0  within  a 
few  meV. 

A  gap  of  similar  size  was  reported  in  a  scanning  tunnel¬ 
ing  spectroscopy  study  of  graphene  [26,27],  and  it  was 
attributed  to  phonon-mediated  inelastic  tunneling  [28]. 
However,  there  have  been  disagreements  about  the  gap 
size  or  even  the  existence  of  gaps  [29-32].  The  inconsis¬ 
tencies  remain  controversial  and  could  be  related  to  sample 
inhomogeneity,  tip  structure,  or  tunneling  at  the  buried 
interface.  The  mechanism  leading  to  the  ARPES  gap  does 
not  apply  directly  to  the  tunneling  case.  However,  for  a  well 
collimated  tunneling  current  perpendicular  to  the  graphene 
layers  from  a  tip  with  a  simple  density  of  states,  the  same 
parity  selection  rules  apply,  and  a  gap  of  —67  meV  can  be 
expected. 

ARPES  is  free  from  tip  effects;  a  further  advantage 
is  angular  resolution.  Experimentally,  the  gap  disappears 
for  increasing  emission  angle  or  in-plane  momentum 
[Fig.  4(a)];  the  symmetrized  spectra  are  for  3  ML  graphene 
on  SiC  at  60  K.  The  two  vertical  lines  indicate  the  edges 
of  the  gap  at  normal  emission.  The  corresponding  first 
derivatives  [Fig.  4(b)]  show  peaks  at  the  edges  of  the 
gap.  The  gap  diminishes  and  gets  filled  in  for  increasing 
emission  angles  because  phonons  at  lower  wave  vectors, 
and  lower  energies,  become  involved. 

While  the  67  meV  phonons  are  the  dominant  contri¬ 
butions  to  the  gap  at  normal  emission,  other  phonons  at 
— 150  me V  are  expected  to  contribute  as  well,  how¬ 
ever  weakly.  The  dipole  transition  term,  even  though 
suppressed,  is  nonzero,  and  its  parity  permits  the  excitation 
of  the  in-plane  phonons.  High-resolution  scans  indeed 
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FIG.  4  (color  online).  Dependence  of  the  gaps  on  k. 
(a)  Symmetrized  spectra  for  3  ML  graphene  as  a  function  of 
in-plane  momentum  k  along  the  T  K  direction.  The  photon 
energy  used  was  50  eV.  (b)  Absolute  value  of  the  first  derivative 
of  (a).  The  zero  level  for  each  spectrum  is  indicated  by  a  tick 
mark  on  the  right  vertical  axis. 


show  a  very  weak  secondary  gap  at  —150  meV  (see  [18]). 
The  results  suggest  the  possibility  of  detailed  phonon 
spectroscopy  based  on  this  method. 

The  present  study  demonstrates  that  ARPES  spectral 
regions  forbidden  for  direct  emission  provide  a  valuable 
test  ground  for  detailed  investigations  of  elementary  inter¬ 
actions.  Electron-phonon  coupling  is  a  weak  interaction, 
but  its  signature  shows  up  clearly  because  of  the  absence 
of  other  spectral  features.  This  insight  is  important  for  a 
detailed  understanding  of  the  basic  physics  of  graphite  and 
graphene  [33,34].  Possible  applications  of  the  method 
include  phonon  spectroscopy  and  determination  of  the 
electron-phonon  coupling  strength  through  a  detailed 
analysis  of  ARPES  data.  Specifically,  the  phonon  gaps  at 
general  points  in  the  Brillouin  zone  could  be  used  to  extract 
detailed  phonon  dispersion  relations,  while  the  spectral 
height  of  each  phonon  gap  provides  a  direct  measure  of 
the  electron-phonon  coupling  strength. 
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Electronic  chirality  near  the  Dirac  point  is  a  key  property  of  graphene  systems,  which  is  revealed  by  the 
spectral  intensity  patterns  as  measured  by  angle-resolved  photoemission  spectroscopy  under  various 
polarization  conditions.  Specifically,  the  strongly  modulated  circular  patterns  for  monolayer  (bilayer) 
graphene  rotate  by  ±90°  (  ±  45°)  in  changing  from  linearly  to  circularly  polarized  light;  these  angles  are 
directly  related  to  the  phases  of  the  wave  functions  and  thus  visually  confirm  the  Berry’s  phase  of  tt  (2tt) 
around  the  Dirac  point.  The  details  are  verified  by  calculations. 
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Electronic  chirality  and  the  Berry’s  phase  are  topics  of 
central  importance  to  graphene  physics  [1-3].  The  direc¬ 
tion  of  the  isospin  vector  in  monolayer  graphene  tracks 
the  rotation  of  the  electronic  crystal  momentum,  resulting 
in  a  Berry’s  phase  of  tt  [4,5].  By  contrast,  the  isospin  in 
bilayer  graphene  rotates  twice  as  fast,  resulting  in  a  Berry’s 
phase  of  2tt  [6,7].  These  unusual  phase  properties  in 
graphene  systems  lead  to  an  unconventional  quantum 
Hall  effect,  which  has  attracted  widespread  interest  in 
utilizing  graphene  as  a  new  platform  for  electronic  appli¬ 
cations  [1].  Angle-resolved  photoemission  spectroscopy 
(ARPES)  has  been  employed  extensively  to  examine  the 
conical  band  dispersion  relations  [8,9].  Previous  studies  of 
monolayer  graphene,  using  either  s  or  p  polarized  light, 
show  that  the  intensity  patterns  have  a  cosine  functional 
form  with  a  maximum  along  the  TK  direction  [9-13]. 
Another  study  found  that  the  intensity  pattern  for  bilayer 
graphene  from  s  polarized  light  has  two  nodes  along  the 
TK  direction,  which  can  be  linked  to  the  Berry’s  phase 
[14]. 

In  this  work,  we  employ  ARPES  with  both  circularly 
and  linearly  polarized  light  to  probe  the  phase  relations  of 
the  wave  functions.  While  these  relations  are  already 
known  from  simple  tight-binding  models,  the  use  of  circu¬ 
larly  polarized  light  allows  direct  sensing  of  the  phases  in  a 
visual  manner.  Specifically,  the  strongly  modulated  circu¬ 
lar  patterns  for  monolayer  (bilayer)  graphene  rotate  by 
±90°  (  ±  45°)  in  changing  from  linearly  to  circularly 
polarized  light.  These  angles  and  directions  of  rotation 
are  determined  by  the  phase  angle  of  the  incident  light 
polarization,  and  arise  from  the  interplay  between  elec¬ 
tronic  chirality  and  coherent  interference  of  dipole  transi¬ 
tions  from  the  s  and  p  polarization  components.  The 
observed  ARPES  patterns  are  in  excellent  agreement 
with  calculations.  The  phases  of  electronic  wave  functions 
are  hard  to  extract  in  general,  but  the  information  is  critical 
for  a  number  of  systems  including,  apart  from  graphenes, 


topological  insulators,  complex  oxides,  and  magnetically 
ordered  materials.  Our  work  on  graphene  layers  as  model 
systems  demonstrates  that  photoemission  with  circularly 
polarized  light  is  a  powerful  approach  to  solve  the  phase 
problem. 

Our  ARPES  measurements  were  performed  at  the 
Synchrotron  Radiation  Center,  University  of  Wisconsin- 
Madison,  using  the  new  U9  VLS-PGM  beam  line.  All 
spectra  were  taken  with  the  sample  temperature  at  60  K. 
The  polarization  purity  is  better  than  99%  for  HP  and  VP, 
and  better  than  80%  for  LCP  and  RCP.  The  accuracy  of 
sample  alignment  with  respect  to  light  polarization  is  better 
than  2°.  Monolayer  and  bilayer  graphene  samples  were 
prepared  on  the  Si-face  of  a  6H  SiC  (0001)  substrate  by 
graphitization  of  the  surface  in  an  atmosphere  of  10“ 6  Torr 
of  disilane  to  ensure  growth  of  graphene  layers  with  large 
domain  sizes  [15,16].  The  coverage  of  graphene  layers, 
excluding  the  6^3  X  6^3  buffer  layer,  is  deduced  from  the 
Dirac-cone  feature  near  the  K  point  [17-19]. 

Our  ARPES  geometry  and  the  coordinate  system  in 
relation  to  the  graphene  atomic  structure  and  the 
Brillouin  zone  are  shown  in  Figs.  l(a)-l(c).  Four  polariza¬ 
tion  configurations  were  employed:  horizontal  polarization 
(HP,  or  p  polarization),  vertical  polarization  (VP,  or  s 
polarization),  left  circular  polarization  (LCP)  and  right 
circular  polarization  (RCP).  The  geometry  is  such  that 
the  vector  potential  A  of  the  incident  light  lies  nearly 
within  the  sample  plane;  thus  Az  can  be  ignored  to  first 
order.  Figure  1  (d)  shows  ARPES  results  around  the  K  point 
measured  along  the  y  direction  with  LCP  and  RCP. 
Evidently,  photons  with  different  helicity  selectively  excite 
one  of  the  two  branches  of  the  Dirac  cone  with  opposite 
isospin  directions.  The  observed  difference  can  be  charac¬ 
terized  by  a  dichroism  coefficient,  D,  defined  as 

D=  ±(/4  — /3)/(/i  +I2),  CD 
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FIG.  1  (color  online).  Experimental  geometry  and  circular 
dichroic  data,  (a)  Schematic  of  experimental  setup.  HP  and  VP 
correspond  to  p  and  s  polarizations,  respectively.  Looking 
toward  the  sample,  LCP  (RCP)  corresponds  to  counterclockwise 
(clockwise)  rotation  of  the  electric  field  vector,  (b)  Coordinate 
system  relative  to  the  atomic  structure  of  a  monolayer  graphene. 

(c)  Surface  Brillouin  zone,  q  is  wave  vector  measured  from  K. 

(d)  Photoemission  spectra  from  monolayer  graphene  near  the  K 
point  measured  along  y  using  30  eV  photons.  The  four  panels 
correspond  to  LCP,  RCP,  their  sum  and  difference. 


where  /  j ,  /2,  /3,  and  /4  are  peak  intensities  of  the  energy 
distribution  curves  as  specified  in  Fig.  1(d),  and  the  +  and 
—  signs  correspond  to  the  7 7*  and  tt  bands,  respectively. 
Measurements  show  that  D  is  close  to  unity  near  the  K 
point  at  30  eV  photon  energy.  This  large  dichroic  effect 
results  from  the  different  characters  of  the  initial  states  on 
the  two  branches  of  the  Dirac  cone;  the  variations  in  the 
photoemission  final  states  are  negligible  over  the  narrow 
range  in  reciprocal  space  [20]. 

Figures  2(a)  and  2(b)  show  typical  constant  energy  maps 
of  the  7 7*  and  tt  bands  of  monolayer  graphene,  respec¬ 
tively,  using  the  four  polarization  configurations.  For  each 
band,  the  semicircular  pattern  rotates  counterclockwise  in 
90°  steps  in  going  from  HP  to  LCP  to  VP  and  to  RCP. 
Furthermore,  the  patterns  of  the  77*  and  tt  bands  are  related 
by  inversion.  The  plots  on  the  right-hand  side  of  Fig.  2 
show  the  measured  intensity  variations  of  the  semicircular 
arc  patterns  about  the  azimuthal  angle  0  of  the  wave  vector 
k  relative  to  the  K  point  [Fig.  1(c)].  The  results  are  well 
described  by  cosine  functional  forms  related  by  90°  off¬ 
sets.  These  angular  offsets  are  related  to  the  phase  of  the 
initial  state  wave  functions,  T^,  which  are  given,  to  first 
order,  by 
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FIG.  2  (color  online).  ARPES  data  from  monolayer  graphene 
taken  with  HP,  VP,  LCP,  and  RCP  using  30  eV  photons.  The 
Dirac  point  is  at  —0.45  eV.  (a)  Constant  energy  maps  at 
—0.05  eV.  (b)  Constant  energy  maps  at  —0.9  eV.  The  right 
panels  show  the  intensity  of  each  semicircular  arc  as  a  function 
of  6.  Each  intensity  curve  has  been  normalized  to  unity  maxi¬ 
mum  intensity. 


'Pf(r)  ^exp(/k  • 

-R<+!*)!  <2) 


where  R;  is  a  lattice  vector,  </>  is  the  carbon  pz  orbital,  a  is 
the  C-C  bond  length,  and  the  summation  is  over  all  lattice 
vectors.  The  +  and  —  signs  correspond  to  the  77*  and  tt 
bands,  respectively.  The  phase  factors  associated  with  the 
two  sublattices,  e±l ^2,  determine  the  chiral  properties  of 
the  electronic  states  and  give  rise  to  a  Berry’s  phase  of  tt 
upon  incrementing  the  azimuthal  angle  6  by  2tt. 

The  allowed  final  state  of  photoemission  near  the  K 
point  must  be  of  even  parity  with  respect  to  reflection  about 
the  vz  plane.  A  straightforward  calculation  yields  the 
photoemission  intensities  for  the  77*  and  tt  bands  [20]: 
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I  oc  | Ax£x(el°/2  ±  e  w /2)  +  Ay£y(el0/2  +  e  l°/2)\2 

=  f  4cos2f  |AX|2|£X|2  +  4sin2§  |Ay|2|£y|2  +  4  sin0Im(Ax£xA*£*), 
1  4sin2f  |AX|2||X|2  +  4cos2 f  lA^ j2 j ^ j2  -  4sin0Im(Ax£xA*£*) 


where  £x  and  gy  are  the  dipole  transition  matrix  elements 
for  the  carbon  pz  orbital  excited  by  the  v  and  y  components 
of  the  vector  potential,  Ax  and  Ay,  respectively.  With 
Ay  =  0  for  HP  and  Ax  =  0  for  VP,  the  first  and  second 
terms  on  the  right-hand  side  of  Eq.  (3)  yield  the  patterns 
excited  by  HP  and  VP,  respectively.  The  results,  propor¬ 
tional  to  1  ±  cos(0),  are  in  good  agreement  with  experi¬ 
ment  [13].  However,  the  results  are  independent  of  the  sign 
of  0.  A  transformation  0^—0  reverses  the  phases  of  the 
chiral  electronic  states  [Eq.  (2)],  but  leaves  the  intensity 
patterns  invariant.  Thus,  measurements  based  on  HP  and 
VP  alone  are  insufficient  to  fully  determine  the  chiral  phase 
relationships. 

This  indetermination  can  be  resolved  by  using  circularly 
or  elliptically  polarized  light,  for  which  the  third  term  on 
the  right-hand  side  of  Eq.  (3)  comes  into  play.  This  term 
contains  a  factor  of  sin(0),  an  odd  function  of  6\  thus,  these 
measurements  should  uniquely  specify  the  sense  of  0.  At 
30  eV,  £x  ~  £y  as  will  be  shown  below.  With  Ax  being  real 
and  Ay  being  imaginary  for  elliptical  polarization,  the 
photoemission  intensity  is  reduced  to 

I  a  2|£x|2Aq[1  ±  cos(0  +  2x)\  (4) 

where  Aq  =  |AX|2  +  |Ay |2  and  x  =  arg(Ax  +  Ay)  is  the 
phase  angle  of  the  light  polarization.  Equation  (4)  shows 
that  the  intensity  patterns  have  the  same  cosine  forms  as 
the  case  of  linear  polarization  except  for  an  angular  offset 
of  2x,  or  twice  the  phase  angle  of  the  elliptical  polariza¬ 
tion.  A  nonzero  offset  removes  the  sign  ambiguity  of  6\  and 
the  patterns  become  asymmetric  with  respect  to  the  trans¬ 
formation  #—►—#.  With  x  =  -tt/4,  0,  77/4  and  77/2  for 
LCP,  HP,  RCP,  and  VP,  respectively,  the  patterns  should 
rotate  clockwise  in  90°  steps  following  this  sequence, 
which  is  indeed  observed  experimentally  (Fig.  2).  The  third 
term  in  Eq.  (3),  resulting  from  interference  between  the 
dipole  transitions  excited  by  Ax  and  Ay,  uniquely  deter¬ 
mines  the  direction  of  rotation. 

The  chirality  of  bilayer  graphene  is  another  interesting 
test  case.  The  Dirac  cone  is  split  into  two  pieces.  For  the 
main  cone  at  the  Dirac  point,  the  phase  factors  e±l6^2  in 
Eq.  (2)  are  replaced  by  e±l°.  Thus,  the  isospin  rotates  twice 
as  fast  as  the  monolayer  case  [7].  As  a  result,  the  intensity 
patterns  become  “compressed”  by  a  factor  of  2  as  a 
function  of  0 .  The  single  intensity  maximum  for  the  mono- 
layer  case  for  each  polarization  configuration  becomes  two 
intensity  maxima  separated  by  180°.  The  angular  offset  for 
the  different  polarization  configurations  should  become 
one  half  of  the  monolayer  case,  or  45°.  Specifically,  the 
two  intensity  maxima  should  be  located  at  0max  =  +x, 
Tx  +  77,  with  the  —  and  +  signs  corresponding  to  the 


I 

7 7*  and  7 7  bands,  respectively.  This  is  indeed  seen  in 
the  experiment  [Fig.  3(a)].  Because  of  the  splitting  of  the 
Dirac  cones,  the  photoemission  intensity  described  by 
Eq.  (4)  should  additionally  be  modified  to  include  this 
effect  [20].  This  explains  the  other  details  of  the  data  and 
why  the  intensity  maximum  at  6  =  77  for  HP  appears 
missing  [14,21].  Figure  3(b)  shows  the  calculated  intensity 
patterns,  which  agree  well  with  the  experiment.  The  char¬ 
acteristic  45°  rotations  of  the  intensity  patterns  between 
linear  and  circular  polarizations,  as  well  as  the  angular 
positions  of  the  intensity  maxima,  are  direct  consequences 
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FIG.  3  (color  online).  ARPES  data  from  bilayer  graphene 
taken  with  HP,  VP,  LCP,  and  RCP  using  30  eV  photons.  The 
Dirac  point  is  at  —0.3  eV.  (a)  Constant  energy  maps  at 
—0.05  eV.  (b)  Calculated  results  with  a  Gaussian  broadening 
of  a  =  0.015  A-1  in  k  space. 
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of  the  electronic  chirality  and  a  Berry’s  phase  of  2tt  in 
bilayer  graphene  (vs.  i r  in  monolayer  graphene). 

The  interference  term  in  Eq.  (3)  depends  on  the  phase 
difference  between  i;x  and  gy.  Let 

Ae,/3  =  €y/gx,  (5) 

where  A  and  (3  are  the  amplitude  ratio  and  phase  differ¬ 
ence,  respectively,  which  can  both  depend  on  the  photon 
energy  because  of  variations  in  the  final  state.  The  ampli¬ 
tude  ratio  A  can  be  extracted  from  the  intensity  ratio 
(=  A2)  between  VP  and  HP  at  6  =  ±tt/ 2  [Fig.  4(a)]. 
The  phase  difference  (3  can  be  extracted  from  the  measured 
dichroism  coefficient  D  [Fig.  4(b)]  via 

The  deduced  (3  [Fig.  4(c)]  is  close  to  zero  for  photon 
energies  above  30  eV,  but  increases  at  lower  photon 


FIG.  4  (color  online).  Deduction  of  the  phase  difference  /3  for 
monolayer  graphene,  (a)  The  photoemission  intensity  ratio 
(=  A2)  between  VP  and  HP  at  ±0.4  eV  relative  to  the  Dirac 
point  as  a  function  of  photon  energy.  The  “error  bars”  indicate 
the  differences  between  results  for  the  77*  and  tt  bands,  (b)  The 
corresponding  dichroism  coefficient,  D ,  obtained  similarly, 
(c)  The  deduced  phase  difference  (/ 3 )  between  transition  matrix 
elements  from  HP  and  VP. 


energies,  where  final-state  band  structure  effects  become 
significant  [20,22].  Our  results  show  that  £x  =  £y  is  well 
satisfied  at  30  eV,  which  leads  to  very  simple  pattern 
rotations  seen  in  Figs.  2  and  3.  At  high  photon  energies, 
band  structure  effects  for  the  photoemission  final  states 
diminish,  and  we  expect  A  — ►  0  and  /3  — ^  0  based  on 
symmetry  considerations  [20].  However,  an  additional 
contribution  from  surface  photoemission  can  come  into 
play  at  very  high  photon  energies  [20,23-25],  and  the 
analysis  can  become  more  involved. 

Our  study  thus  shows  that  polarization-dependent 
ARPES  can  be  a  powerful  tool  for  probing  the  phases  of 
electronic  wave  functions  in  solids.  While  circular  dichroic 
effects  are  well  known  in  gas  phase  ARPES  [26-28],  much 
less  has  been  studied  for  solid-state  systems  because  of 
complications  including  electron  correlation  [29,30].  The 
Dirac  Fermions  in  graphene  systems  are  of  special  interest 
because  of  their  chirality  and  nontrivial  Berry’s  phases, 
which  can  be  conveniently  probed  by  elliptically  polarized 
light  with  a  built-in  rotation  of  the  vector  field.  Interference 
of  the  dipole  transitions  gives  rise  to  the  observed  pattern 
rotations  governed  by  the  ellipticity  of  the  incident  light. 
This  understanding  and  the  methodology  presented  here 
should  be  useful  for  exploring  other  systems  and  problems 
including  topological  insulators,  spin  splittings  and 
magnetic  ordering.  In  particular,  the  surface  states  on 
topological  insulators  should  have  a  nontrivial  Berry’s 
phase,  and  a  recent  photoemission  study  indeed  revealed 
interesting  circular  dichroic  effects  [31,32].  Further  analy¬ 
sis  of  the  topological  states  based  on  the  present  method  is 
underway. 
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